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Abstract 

Although there are phenomenological indications that the low-energy con- 
stants in the chiral lagrangian may be understood in terms of a finite number of 
hadronic resonances, it remains unclear how this follows from QCD. One of the 
arguments usually given is that low-energy constants are associated with chiral 
symmetry breaking, while QCD perturbation theory suggests that at high en- 
ergy chiral symmetry is unbroken, so that only low-lying resonances contribute 
to the low-energy constants. We revisit this argument in the limit of large N c , 
discussing its validity in particular for the low-energy constant L$, and conclude 
that QCD may be more subtle that what this argument suggests. We illustrate 
our considerations in a simple Regge-like model which also applies at finite N c . 
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1 Introduction 



The low-energy physics of the Goldstone bosons associated with the spontaneous break- 
ing of chiral symmetry is described by the chiral lagrangian, the construction of which is 
solely based on the general principles of quantum field theory and the symmetries of the 
underlying fundamental theory of the strong interactions, QCD. The only "dynamical" 
information present in the chiral lagrangian connecting it to QCD is encompassed by 
the values of the low-energy constants (LECs) that parameterize the terms appearing 
in the chiral lagrangian. At lowest order, there are two LECs, the pion decay constant 
f n and the quark condensate (ipip), whereas at higher orders, a growing number of 
LECs appear. In the theory with three flavors, there are ten physical LECs at order 
p 4 [T], denoted as Lj, % — 1, . . . , 10. 

It is thus an interesting question to ask what physics of the strong interaction leads 
to the values of the LECs observed in nature. Since the physical spectrum contains 
infinite towers of resonances with given quantum numbers (spin, flavor and parity), it 
is natural to think of the values of the order p 4 (and higher-order) LECs as resulting 
from integrating out these resonances J2JI3]- 1 In particular, sum rules have been derived 
relating the order-p 4 constants to the masses and decay constants of these resonances. 

Here we revisit in particular the sum rule relating L§ to the scalar and pseudoscalar 
resonances with isospin 1, working in the large- N c and chiral limits. The choice of 
taking N c large simplifies the discussion because hadronic resonances are stable in 
that limit. Furthermore, it is natural to work in the chiral limit, since the LECs are 
independent of the quark masses. For comparison, we will also discuss the sum rule 
relating Liq to vector and axial- vector resonances with isospin 1. 

We begin with reviewing the usual type of sum rule relating L$ to resonance param- 
eters through the difference of scalar and pseudoscalar two-point functions IIs(g 2 ) — 
Ilp(g 2 ), and show that it is in principle ill-defined. We show that the problem originates 
from the fact that each of these two-point functions is quadratically divergent in QCD. 
We will then derive a new sum rule, starting from a subtracted dispersion relation, in 
Sec. 2.2. In Sec. 3 we give an example of how the new sum rule works, using a simple 
Regge-like model in which the exact answer is known. This model allows us also to 
consider the situation at finite N c . In Sec. 4 we test the phenomenological approach in 
which the difference n,s(g 2 ) — ITp(g 2 ) is saturated by a finite number of resonances. The 
usual argument for this approach is that this difference vanishes in perturbation theory, 
suggesting that the contribution from resonances higher in the spectrum cancels, at 
least approximately. Applying this approach to our model, we find that the effective 
parameters (i.e. masses and coupling constants) which are required to fit high-energy 
constraints as well as to yield the correct value for L 8 are not the same as those of the 
original infinite spectrum. The final section contains our conclusions. 

lr The physics of and (ipip) is associated with the breaking of chiral symmetry by the vacuum. 
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2 The large- A^ c limit 

2.1 Review 

In the large- N c limit, L 8 is given by @] 

i6B 2 L 8 = An(o) = l - (n s (o) - n'p(o)) , (1) 

in which the chiral limit is taken, and where we define 

n 5 ,p(<7 2 ) = i f d 4 x e*** (0\t{JsAx)4A°)} l°> > ( 2 ) 

with Js(x) = d(x)u(x) and Jp{x) = d{x)i^u(x). The order-p 2 constant B is equal to 
— (ipifj) I / 2 in the chiral limit. The prime on ITp(O) in Eq. (JJJ indicates that we omit the 
pion pole present in the pseudoscalar two-point function Up(q 2 ). Since the difference 
Us — lip cancels in perturbation theory for vanishing quark masses, this means that 
AII(g 2 ) obeys an unsubtracted dispersion relation in the chiral limit, i.e. 

Antf > 3 I (iW) - n^)) = ^ f * iSLg^zaW) . (3) 

Z Z7T ^/q 6 — y — 6c 

Equation (JTJ) states, then, that L 8 becomes a physical quantity at iV c = oo. In this 
limit, there is an infinite tower of infinitely narrow resonances in each channel, and we 
thus have that |5] 

-y OO 

-Imris(i) = 2y>J(n)<y(*-Aff(7i)), 

7T z — ' 

n 

1 oo 

- Im U' P (t) = 2 V Fp(n) 5(t - Mp(n)) . (4) 

7T f ^ 



Use of Eq. (J3J) in Eq. (JHJ) leads formally to 

Allfa 2 ) = V Fg2(n) V F ^ (n) (5) 

[q ] "' ^ M 2 (n) - g 2 - «e ^ M 2 (n) - q 2 - le ' 1 ' 

In arriving at this expression for AII(g 2 ), we assumed that the integral over t and 
the sum over n may be freely interchanged. Ignoring the question of validity of this 
assumption for the time being, one thus arrives at a sum rule for L§: 



As it stands, this sum rule is not well-defined, because both sums appearing on the 
right-hand side are quadratically divergent. First, let us give an example. If we assume 
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Regge-like behavior asymptotically for large n, with 2 



M s,p( n ) ~ nA > n lar g e > ( 7 ) 
F| P (n) ~ «:A 2 M| P (n) , 

with A a constant of order 1 GeV, and k is equal to the coefficient of the parton model 
logarithm, then Eq. © looks like 

(oo oo \ 

n n / 

which is clearly ill-defined. However, as we will now discuss, the problem is model 
independent. We will return to the model (J7J) in Sec. 3. 

In order to proceed in the general case, let us introduce cutoffs N$ t p on the sum- 
mation index n in the scalar and pseudoscalar channels, corresponding to cutoffs 
Ms(Ns), Mp(Np) on the masses, in order to regulate the sums. Furthermore, we 
require that both two-point functions match on to QCD perturbation theory for asymp- 
totically large euclidean Q 2 = —q 2 . In particular, the leading large-Q 2 behavior is given 
by nQ 2 logQ 2 , with k a constant common to both channels, because the perturbative 
expressions for scalar and pseudoscalar two-point functions, valid for asymptotically 
large Q 2 , are equal. It follows that indeed the number of resonances in each channel 
has to be infinite j^j. Furthermore, writing 

F 2 s, P (n) = fl P {n)Ml P {n) , (9) 

each of the sums in Eq. (jSJ) can be written as (omitting the indices S or P) 

Y F ^ n) = Tf 2 (n)-o 2 T p[n) (10) 



f\n) 



dM 2J v 7 M 2 + Q 2 ' 



If (dn/dM 2 )f 2 (n) goes like a constant for large M 2 , the second integral will behave like 
log (M 2 (N)/Q 2 ) for M 2 (N) S> Q 2 , thus yielding the desired parton-model logarithm. 
As a consequence, the first integral in Eq. (JTUJl is quadratically divergent. Note that 
M 2 (N) ^> Q 2 should hold for any Q 2 , and that we thus need to require that M 2 (N) — > 
oo when iV — > oo. It follows that the constant k can be expressed as 

K = 7^2 fs( n ) = lim -TTH J pin) ■ (II) 

n-+oa dMg n—>co amp 

Furthermore, it also follows that Mg(Ns) / 'Mp(Np) — > 1, to avoid a term linear in Q 2 
mU s (q 2 )-U P (q 2 )Jorq 2 = -Q 2 . ' 

2 For the asymptotic proportionality between decay constants and the masses, see Eq. (Illf) below. 
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In fact, we may choose to introduce a common cutoff A 2 Q on the t integrals in Eq. (J3J), 
which implies cutoffs Ng t p on the sums over n such that Mg(Ns), Mp(Np) < A 2 D , while 
Mg(N s + 1), Mp(N P + 1) > A 2 D . Using Eq. (0), our regulated sum rule now reads 

N s (A co ) N P {A co ) 

16B 2 L 8 (A co ) = fs(n)~ £ f P {n). (12) 

n n 

This is of course finite, by construction. However, it does not satisfy the basic field- 
theoretic requirement that physics should be independent of the cutoff, i.e. that the 
relative change in L 8 due to a change in cutoff goes to zero as we take the cutoff to 
infinity. Mathematically, this requirement translates into 

L s {A co )=L B + o(J-\ , (13) 

where L 8 is the value of the LEC in the limit of infinite cutoff. That Eq. ()12|) does 
not in general satisfy this requirement can be seen as follows. When we increase A co , 
both Ns and Np increase by integer steps when A co moves past the next resonance 
mass in either channel. Suppose that we increase A co such that Ns changes by 1 to 
N' s = N,s + 1, while N' P = Np stays the same. This leads to a change AL§ in L 8 equal 
to 

Ai8 .fc^±i). (14) 

We do not at present have enough information about QCD allowing us to conclude 
that the right-hand side of this equation goes to zero for A co — > oo. In fact, it is not 
the case for the Regge-like behavior of Eq. ((Zj), for which /|(n) goes to a constant for 
large n (cf. Eq. (|TT|) ). Therefore, we conclude that the sum rule for L 8 as given by 
Eq. ©, even with the cutoff we introduced to regulate the sums in Eq. (0), is not 
well-defined. 3 



2.2 A better sum rule for L§ 

Obviously, the invalid step in the derivation of the sum rule © is the interchange of the 
integral over t and the sum over n, as this led to an expression containing the difference 
of two quadratically divergent sums. We may remedy this problem by starting from a 
once-subtracted dispersion relation for AII(g 2 ) [Uj, 

Q2 roc llm(Us(t)-U' p (t)) 

auV = -q 2 ) = An(o) dt ; ; (15) 

where again the pion pole has been omitted. Substituting the large-iV c expressions 
Eq. (0J) for the spectral functions, one obtains 

Antf = -tf) = An,o, - Q > (± - £ J £fU , (i.) 



'Unless there is an N so that fs(n) = fp(n) for all n > N. 
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where we used Eq. (JSJ), and where Ns and Np are to be taken to infinity. Again we 
interchanged the integral over t with the sum over n, but now each sum in Eq. (|16j) is 
only logarithmically divergent: 

t ~ dM2 (^ /2(n) ) WTc? ~ K log {mQ2) • (17) 

(cf. Eq. (JUJ)), and therefore a cutoff A co needs to be introduced such that Q 2 -C A 2 D ~ 
Mf(iV s ) ~ Mp(Np). The limit of N s and iV P to infinity thus has to be taken in 
a correlated way such that this condition is fulfilled. It follows that the difference 
in Eq. (jTSj) is finite and unambiguous. Using the fact that AU(Q 2 ) should vanish for 
asymptotically large Q 2 , we arrive at a fully regulated sum rule for L 8 , valid in large- N c 
QCD: 4 

16B 2 L 8 = AII(O) = lim Q 2 lim ( V ^ 2 - V A J?[ n) —] , (18) 

v ' q 2 -*oo 7v s ,iVp^oo M|(n) + Q 2 z ^M 2 (n) + Q 2 J v ; 

where it is understood that the limit iVs, Np — > oo is taken in a correlated way, such 
that M 2 (N S )/M 2 (N P ) -> 1. 

If, as already discussed in Sec. 2.1, we take the number of resonances in each channel 
to be finite, assuming that there are ng in the scalar channel and rip in the pseudoscalar 
channel below a certain scale sq, above which we stipulate that the scalar and pseu- 
doscalar two-point spectral functions exactly cancel, each of the sums in Eq. (J18)) is 
finite, and we may take the large-Q 2 limit under the sums, obtaining 

n S n P 

16B 2 L 8 = ^f 2 s (n)-J2fUn) , (19) 

n=l n=l 

i.e. a finite version of Eq. (JHJ). In general, if the difference of the sums in Eq. fll8j) would 
be sufficiently convergent, one would be allowed to take the large-Q 2 limit under the 
sums. However, it is known that the coefficients of the OPE in QCD have anomalous 
dimensions, implying that logarithmic corrections to the inverse powers of Q 2 appear. 
This already implies that expanding in 1/Q 2 under the sums is in general not allowed. 
In particular, an exact cancellation of the spectral functions above some So cannot 
occur. 

It is instructive to compare our sum rule for Lg with a similar sum rule for Lxo, 
which is related to the difference of the vector and axial-vector two-point functions 
Uv{q 2 ) and ILa (<? 2 ) at q 2 = 0. The key difference is that in this case each of these 
two-point functions is only logarithmically divergent to begin with, because of gauge 
invariance. The vector and axial-vector two-point functions are defined by 

U^ A = tj dx (0\T { jrjx) jj^(0)} |0> = (q»q„ - g, u q 2 )U VA (q 2 ) , (20) 
The pion pole does not play any role in this sum rule. 
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with Jy(x) = d(x)j II u(x) and J A (x) = d(x) r y fl r y 5 u(x). Perturbation theory tells us that 
these two-point functions should both behave like dlogQ 2 for asymptotically large Q 2 , 
with k' another constant fixed by the parton model logarithm for this case. For large 
N c , the spectra look like 

1 oo 

-lmU v (t) = 2V F v (n) 8(t - M v {n)) , (21) 

n 

1 oo 

-Imn^(t) = 2^Fl(n)^-M»), 

n 

where again the prime indicates that we omit the pion pole. A line of reasoning similar 
to that of Sec. 2.1 shows that 



dn „ , s , dn 



K ' = lim JMl F v {n) = lim -^ F ^ n) ■ (22) 

ra— >oo ClMtr n— »oo CLM . 



From this it follows that both II y and 11^ are logarithmically divergent, so that no 
subtraction analogous to that in Eq. ([T5]) is needed. We thus obtain the sum rule 

-4L 10 = i(n v (o)-n^(o)) (23) 



2 



n v ,n a ^oo Af£(n) ^ A^IH 

where the limit Ny,N A — > oo is again taken in a correlated way such that i\ 2 co ~ 
M V (N V ) ~ M_|(JVa) I3IH1- Now, if we increase A co such that, say, N v ^ N v = N v + 1 
while N A ->N' A = N A , L 10 changes by AL W = -F v (n = N v + l)/AM v (n = N v + 1). 
If we assume that M v (n) grows like n (as in Regge-like behavior), we have that ALio 
1/Ny — > for Ny — > oo. This is to be compared with the case of L§ in Eq. ()14j) . 



3 A Regge-like model at finite N c 

Clearly, the infinite sums in Eq. (@J) are not well-defined mathematical expressions. 
However, at finite N c the Dirac delta distributions should become better behaved, and 
it is therefore illustrative to consider the case of finite N c . In order to do this, we 
will consider a simple Regge-like model. In this model, each scalar resonance will be 
represented by a pole of the form [HI ED] 

2F 2 (n) , x 

(24) 



zA 2 + M 2 {n) ' 
with 
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and similar for the pseudoscalar resonances, with masses and residues given by 



M z SP {n) = m 2 SP + nA 2 , n = 0, 1, . . . , (26) 



F 2 (n) = nk 2 M 2 (n) 



with ms,p also of order 1 GeV (cf. Eq. (|7jl ). For simplicity, we take the parameters A 
and a the same in both channels. This choice satisfies Eq. ([lip. 

Expanding the denominator of Eq. (|24|) to leading order in 1/N C , there are poles 
near q 2 = M| p (n)(l — ia/N c ). Our model thus describes resonances with masses 
~ M s> p{n) and decay widths 



^s,p 



,(n) ~ aM SiP (n)/N c . (27) 



In the limit N c — > oo, all resonances are stable. The function (J24|) has the correct 
analytic behavior in the complex q 2 plane jU] : it is analytic everywhere on the physical 
sheet, except for a cut along the positive real axis, starting at q 2 = [9 . This model 
does not include all corrections one expects at finite N c ; in particular, it does not 
include the multi-particle continuum starting at q 2 = due to coupling of the sources 
Js,p to pions, which would make L 8 run. However, as we will see below, the included 
finite- N c behavior already leads to an interesting observation. 

The scalar spectral function for our model is given by (for simplicity, we set the 
scale A equal to one in the rest of this section) 5 

*ti(t< cos(Cir) + M») 2 + (f< sm(Cir) + e) 2 

which for A" c — ► oo, i.e. £ — > 1, reproduces Eq. (J3J. The expression for Im ITp(t) is 
similar. These spectral functions are logarithmically divergent with n, and we thus 
introduce a cutoff N = iVs = ATp + c with c a finite arbitrary constant, so that 
Mg(Ns) I Mp(Np) — > 1 for — > oo. 6 The difference between the scalar and pseu- 
doscalar spectral functions is well defined, and we have that, starting from the disper- 
sion relation in Eq. (jSJ), 

An«» = = r~t( ( " + m|)it T (C "' - 2 -(s^P)) , (29) 

Wo t ^ ^(«<cos(Cx)+n + m|) 2 +(*tsin(C7r)) ! / 

where we set e = 0, which we are allowed to do as long as we keep A^ c finite. We may 
also send N c to infinity, but in that case we need to keep an infinitesimal e in order to 
reproduce the Dirac 5-functions of Eq. (j3J). 



5 Note that we start our sums at n — 0, whereas in Sec. 3 of Ref. sums start at n — 1. 
6 We may choose N$ = Np + c with c an arbitrary finite constant because nothing will depend on 
c in the limit of infinite cutoff 0. 
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Substituting y = tt and performing the sum over n, we obtain 



16B 2 L 8 



An(o) 




/ dy (e- ln ^( ye - 1 ^ + m |) _ e ^^y^ + m |) _ (5 



-r oo 



J 



(30) 



where 



d 
dz 



logT(z) 



(31) 



is the digamma function. 7 Interestingly, in this simple model, the result is independent 
of N c , and it is thus also valid for N c — > oo. 8 The technical reason for this is that 
setting Q 2 = in order to obtain L 8 removes the dependence on N c , which only 
appears through the exponent (. We note that the new sum rule Eq. (|T5j) does apply 
to our model at finite N c as well, if one simply replaces Q 2 by zA 2 with z as defined in 
Eq. ()25|) . and then takes the limit z — > oo in Eq. (j!8|) . 

Let us comment on this result. First, if we calculate L 8 for this simple model using 
our new sum rule, Eq. (|18|) . we obtain exactly the same result. 9 Second, a key point 
in the calculation of Eq. (j3*Uj) is that we first performed the sum over n, and after that 
the integral over t. These operations do not commute: if we would interchange the 
sum and the integral, the result would be different. If we substitute y — in Eq. (}2Tjj) 
and integrate term by term over y, we obtain for each term 



Clearly, the result (|3^j) does not equal the correct result, Eq. (pJj) . Moreover, it depends 
on the undetermined constant c, i.e. on the details of the regulator. If, as in Sec. 2.1, we 
increase A co , c alternates between the values and 1 (taking m| < m 2 p ), and the value 
of 16-B 2 AL 8 defined in Eq. (I14|) between kA 2 and 0. This simple example calculation 
demonstrates precisely what goes wrong in the derivation of the naive sum rule Eq. (jHJ): 
as we already noted in Sec. 2.1, that equation is obtained by first performing the t 
integrals in Eq. (jSJ) (without subtraction), and then the sums over n, thus leading to 
the difference of two quadratically divergent sums. The derivation of the new sum rule 
(I18|) avoids this problem by starting from a subtracted dispersion relation, Eq. (|15|) . 

7 Our explicit calculation refutes the claim of Ref. ^] that Lg vanishes in this type of model. 
8 We would obtain the same result setting N c = oo from the outset, but taking e — > at the end of 
the calculation. 

9 This calculation was done in Sec. 4 of Ref. [H]- 




(32) 



and, restoring A, one would thus find that 




(33) 
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Figure 1: Plot of y times the function in Eq. \SU\) with m 2 s = (0.8167) 2 , mp = (1.083) 2 and 
( = 0.82, as a function of y in units in which A = 1, with arbitrary units on the vertical axis 
(solid line). The dashed line is the OPE to leading order. 

Finally, we observe that in our model the result is independent of N c , and thus it is 
equally valid for infinite as well as finite N c . In this respect, it is interesting to consider 
the spectral function at finite iV c , which we plotted in Fig. 1. More precisely, Fig. 1 is 
a graph of y times the integrand of Eq. (JBUj) as a function of y, for a value of ( close 
to unity. We see that the asymptotic behavior predicted by the OPE sets in already 
before y = 3, and that the resonances at higher values are completely "washed out" 
by their growing decay widths. In contrast, at infinite N c the same plot would show 
an infinite sequence of Dirac delta functions, cf. Eq. (HJ). Nevertheless, despite this 
very big difference in the behavior of the spectral functions at finite and infinite N c , 
the value of L 8 is the same in both cases. Even though the perturbative tail cancels 
out in the difference between the scalar and pseudoscalar channels, L 8 is not saturated 
by a finite number of resonances in this model. The contribution to the integral in 
Eq. (f3T)j) coming from the region where the asymptotic behavior already sets in is 
equally important as that from the region where resonances are clearly visible. 

Our example model can be made more realistic by including a finite number of 
resonances in the scalar and pseudoscalar channels below a certain scale so which are 
not part of the equally-spaced towers of Eq. pij) (so that m| P > so). In that case, 
with ris(p) extra states in the scalar (pseudoscalar) channel, we would obtain 

n s np 

16B 2 L 8 = ]T f$(n) - fp(n) + <rn 2 P - m 2 s ) . (34) 

n=l n=l 

This is to be compared to Eq. (|19j). The extra term in Eq. originates (in our 
model) from the infinite towers present above sq. While we do not know what this 
extra term will look like in the case of QCD, we do know that at large iV c there is an 
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infinite number of resonances in each channel. Therefore, the omission of a term like 
this in Eq. (|19|) would represent an unknown systematic error in the application of this 
sum rule to QCD. 

While we certainly do not expect our model to be very realistic, it turns out to be 
interesting to apply it to a numerical example, already in the large- N c limit. In order 
to do this, we first want to eliminate the unknown ratio k/B 2 from Eq. (j3*Uj) . Requiring 
absence of a 1/Q 2 term in the OPE expression for AU(Q 2 ) in this model, one obtains 
the additional relation [U| 10 

2B 2 f 2 = K{m 2 s + m 2 P - A 2 ) (m 2 P - m|) , (35) 

where / = f n in the chiral limit. Combining this with Eq. (J3U|) . we find 

r 1 P 



m 2 s + m 2 p — A 2 



(36) 



Taking m s = 0.98 GeV, m P = 1.3 GeV, A = 1.2 GeV (determined from the a (1450) 
and 7r(1800)) and / = 87 MeV, we find L$ = 8 x 10~ 4 to be compared with, e.g., 
L 8 (M P ) = 9 x 10" 4 [IE]. 

Let us again contrast the case of L§ with that of Lio, as we already did at the end 
of Sec. 2. In a model with an equally-spaced tower of resonances in each channel, one 
finds for L w I7j 



-AT = ;r'V^ \ m A-" b V)^ ( o 7 \ 

^ (™v + nA 2 )(m 2 A + nA 2 ) ' { } 

in which my,A are the lowest masses in the two towers. The key difference is that in this 
case the contribution from the towers is not only proportional to the difference m 2 A —my, 
analogous to the case of L 8 , but also suppressed by inverse powers of the product of 
these resonance masses. Taking the contribution from the first few resonances in the 
tower may already be a good approximation. For instance, with k' = l/8n 2 and 
taking m v = 0.77 GeV, m A = 1.25 GeV and A = 1.3 GeV gives L 10 = -5.6 x 10~ 3 , 
to be compared with the full sum in Eq. ()37|) which gives L w = —7.1 x 10~ 3 . The 
experimental estimate is L 10 (M p ) = —5.5 x 10~ 3 [Tfij . 



2 - m^A 2 



4 Saturation with a finite number of resonances 

It is often assumed that each spectral function may be approximated by a finite sum 
over resonances up to a certain scale so, and a "perturbative continuum" at values 
of t > s . The perturbative continuum then cancels between scalar and pseudoscalar 
channels, leaving the difference between two finite sums in Eq. (|T^|) . In these sums, the 
decay constants and masses of the resonances are to be considered effective parameters 
to be fixed by requiring consistency with the low- and high-momentum expansion 
of the relevant Green's function, as obtained from the chiral and operator product 

10 Here ms,p correspond to the lowest masses in the towers, whereas in Ref. their values were 
expressed as Jm 2 s P + A 2 . 
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expansions, respectively In the large- N c limit, this phenomenological ansatz is a 
rational approximation j^J to QCD Green's functions, known as the Minimal Hadronic 
Approximation (MHA) [3 . It encompasses vector meson dominance as a particular 
case. 

One may also attempt to use only the high-energy input from the OPE to predict 
the low-energy constants, which parameterize the low-momentum expansion. Although 
phenomenologically this assumption seems to work it remains unclear how this 
follows from QCD. It is also important to be able to control the systematic error 
associated with this approximation since, in practice, it is very difficult to go beyond 
the inclusion of one resonance per channel. Furthermore, there are often not enough 
equations to determine the decay constants and the masses, so that the effective masses 
are actually guessed from the position of the physical masses. 

In order to gain some insight into these questions a model like that in Eqs. 



is useful because it allows for a comparison of the phenomenological ansatz with the 
exact large- iV c results produced by the model. Therefore let us consider the ansatz 

l< n * !) - n ' ( ' ! »-=?fsj-^-?' (38) 

where Ms,p and Fs,p are parameters to be determined. The idea is to use short-distance 
(OPE) constraints to fix these parameters, and then use their values to predict Lg from 

16B 2 L S = (39) 



M| M 2 



Writing 



Fl,p = flpMl P , (40) 
the leading OPE {l/Q 2 ) constraint is [S] 

/|M| - f 2 P M 2 = ~(r4 - m|)(m| + m 2 P - A 2 ) . (41) 

With values of m^p and A close enough to 1 GeV, and assuming that nip > nis, as 
phenomenology suggests, the right-hand side of this equation is positive. In fact, this 
is what happens for the values chosen after Eq. (}3T)j) . 

While this one equation is not sufficient to determine the MHA parameters, an 
interesting observation can already be made: MHA does not coincide with the physical 
masses and decay constants. Suppose that we look for solutions in which the decay 
constants take their actual values. In our model, that means /| = fp — k,A 2 . One 
immediately concludes that the ordering of the mass parameters comes out reversed, 
i.e. Ms > Mp while the true ordering is mp > ms- The mathematical origin of this 
difference in sign comes from the fact that the model contains an infinite number of 
resonances. 

To investigate MHA in more detail, let us impose that it gives us the right value of 
L 8 , which adds the relation 

fl-f 2 P = K(m P -m s ) . (42) 
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The lowest resonance masses are relatively well known, 11 while the / parameters are 
much less well known. If one takes the MHA masses as the physical ones, i.e. M$ = ms 
and Mp = mp, one can solve for f$ and fp to find 

fl = ^(m P -m| + A 2 ) , (43) 
fp = ^{m 2 s -m 2 P + A 2 ) . 

These solutions are very far from their "real-world" values (which, in our test, are the 
model values fs = fp = k-A 2 ). Taking the ratio, and using for instance the numerical 
example discussed at the end of Sec. 3, we find 

fl ml-ml + K 2 ^ 1 

fl m 2 P - m 2 s + A 2 " 3 ' 1 ! 

instead of 1, which is the actual value. Note that both the low-energy constraint (}4*2*|) 
and the high-energy constraint (}4~Tj) force /J to be different from fp to the extent 
that ms is not equal to mp. Remarkably, nevertheless, the MHA expression (pIHj) 
approximates the true function Ilg — lip to within a few percent for all euclidean Q 2 . 
Similar properties of the Adler function were studied in more detail in Ref. [T2] . 

One may also try an even simpler ansatz, in which only one scalar resonance (in 
addition to the pion) is kept. Solving Eqs. (j41j) and ()42|) for /J and Mj, we find, using 
the same numerical example, M| ~ 0.63 m| and /| ~ 0.51 kA 2 . As before, the MHA 
values are very different from the real ones, even though the rational approximation to 

— n P works again within a few percent for all euclidean Q 2 . 



5 Conclusion 

In this paper, we revisited the connection of the order-p 4 LEC L 8 to the scalar and 
pseudoscalar resonance parameters. In Sec. 2.1 we reviewed the fact that a sum rule of 
the form © is generally ill-defined, and showed that this is to be expected on general 
grounds in QCD. In Sec. 2.2 we derived a better sum rule, Eq. (fTSj) . and argued why 
this new sum rule gives a finite and universal result. Universality follows because 
each of the sums in Eq. (|18|) is only logaritmically divergent, with the difference being 
finite, and thus the precise details of the choice of the cutoffs Np relative to Ns do not 
matter, as required by quantum field theory [8 j . For instance, increasing Ns to Ns + 1 
in Eq. ([!%]) does not change L$ in the limit Ns — > oo. 

Of course, the ill-defined sum rule of Eq. © itself is never used in practice, but 
instead each sum is restricted to a relatively small set of resonances below a scale s 
with adjustable parameters, as reviewed in Sec. 4. The derivation of such sum rules, 
in this case Eq. (|19p. assumes that above a certain scale so the spectral functions 
Im Hs(t)/7r and Im IIp(i)/7r are exactly equal. In contrast, our new sum rule, Eq. (fTSj) 

11 Assuming that the lowest scalar surviving the large- N c limit is the ao with mass 0.98 GeV. 
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relies only on the fact that II5 and lip become equal for asymptotically large Q 2 in 
the euclidean regime. 

In Ref. [Tl] it was shown that, by using the theory of Pade approximants there 
is a mathematically well-defined connection between resonance parameters and LECs 
in large-iV c QCD. However, this work is not applicable in the case of L§ because the 
spectral function defining Afl(g 2 ) is not positive. 

The sum rule (II 9j) can approximate the sum rule (|18|) if the contribution from the 
resonances above a certain sq is very small. However, there appears to be no reason 
for this to be true in QCD, unless some new mechanism plays a role higher up in the 
spectrum. A possibility might be that the resonance parameters of S and P states 
pair up at higher scales, aligning sufficiently fast for higher resonances to almost not 
contribute to our sum rule, Eq. (fTHjl. This scenario of so-called "parity doubling" or 
"chiral symmetry restoration" has been the subject of much speculation recently, but 
little to nothing is known to date about its validity |15j . In fact, our simple model 
demonstrates that this parity doubling in QCD does not follow from the cancellation 
of the perturbative tail in Us — Up, and neither does the saturation of L$ by a finite 
number of resonances. 

Clearly, our new sum rule is less practical if indeed the contribution from resonances 
above sq cannot be ignored. However, ignoring this contribution anyway would obvi- 
ously introduce an unknown systematic error. In Sec. 3 we showed how the new sum 
rule works in a simple example, in which an explicit form for the contribution coming 
from resonances above so can be derived, cf. Eq. (JHH). While Sec. 2, and in particular 
the sum rule of Eq. (jl8)l . apply to QCD in the limit of infinite N c , our model is defined 
for any value of N c . It is semi-realistic in the sense that it has consistent analytic 
behavior in the whole complex q 2 plane, while it predicts widths which are suppressed 
like 1/N C , but which grow with increasing resonance mass. In particular, we point to 
Fig. 1, which shows a qualitatively realistic spectral function for the two-point function 
An(g 2 ). Interestingly, the result for L 8 in that model, Eq. (J3U|) . is independent of N c . 
We thus speculate that also at finite N c our general conclusion, that the region above 
So cannot be ignored, remains valid. 

Phenomenologically, Eq. (}3*9"j) seems to work quite well [T^]- As we have argued in 
Sec. 4 using our Regge-like model, this may be at the expense of values for the pa- 
rameters in the rational approximant (J38|) . which substantially deviate from the actual 
ones representing the lowest resonances, but which produce a very good approxima- 
tion to the function lis — Hp over the entire euclidean region in Q 2 . Since we lack 
precise information on the scalar and pseudoscalar decay constants in the real world, 
such deviations might go unnoticed. We conclude that, in spite of the fact that the 
phenomenology of the scalar/pseudoscalar sector looks very reasonable, a deeper un- 
derstanding of the relation between L 8 and the spectrum in the scalar and pseudoscalar 
sectors in large- N c QCD remains an interesting puzzle. 
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